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Abstract. We present a special science case for an intermediate resolution multi- 
^ ■ object IR spectrograph for the VLT. We have constructed new models for massive 

Population III stars and metal-free stellar populations (see Schaerer 2001). The prop- 
erties of individual stars and integrated populations, including their ionising fluxes and 
SEDs are discussed. We also study their dependence on the poorly known IMF and 
the star formation history. The synthetic spectra are used to simulate spectroscopic 
observations of the expected emission lines of Pop III galaxies. We show that such an 
^ ' instrument should allow the detection and efficient observations of the first galaxies 

Tj" ' with the VLT. 
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Important advances have been made in recent years on the modeling of the first 
stars and galaxies forming out of pristine matter — so called Pop III objects 
- in the early Universe (cf. review of Loeb & Barkana 2001, proceedings of 
Weiss et al. 2000 and Umemura & Susa 2001). With no doubt these efforts are 
motivated by the approaching possibility of direct observations of such objects 
at very high redshift with NGST and large ground-based telescopes. 
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^ What are the expected observational signatures of Pop III galaxies ? Tumlin- 

son & Shull (2000, hereafter TS00) have recently pointed out the exceptionally 

strong He+ ionising flux of massive [M > 40 M Q ) Pop III stars, which must be a 

/\i • natural consequence of their compactness, i.e. high effective temperatures, and 

non-LTE effects in their atmospheres increasing the flux in the ionising continua. 
As a consequence strong He n recombination lines such as He n A1640 or He n 
A4686 are expected; together with AGN a rather unique feature of metal-free 
stellar populations, as discussed by TS00 and Tumlinson et al. (2001, hereafter 
TGS01). Instead of assuming a "standard" Population I like Salpeter IMF up to 
100 M Q as TS00, Bromm ct al. (2001, hereafter BKL01) have considered stars 
with masses larger than 300 M Q , which may form according to some recent hy- 
drodynamical models. An even stronger ionising flux and stronger H and He n 
emission lines were found. 

Some strong simplifying assumptions are, however, made in the calculations 
of TS00, TGS01, and BKL01. 
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1) All stars are assumed to be on the zero age main sequence, i.e. stellar evolu- 
tion proceeding generally to cooler temperatures is neglected (but cf. TGS01 
for a simple estimate of this effect). 

2) BKL01 take the hardness of the stellar spectrum of the hottest (1000 M©) 
star as representative for all stars with masses down to 300 M©. This leads 
in particular to an overestimate of the He n recombination line luminosities. 

3) None of the studies include nebular continuous emission, which cannot be 
neglected for metal-poor objects with such strong ionising fluxes. This pro- 
cess increases significantly the total continuum flux at wavelengths redward 
of Lyman-a and leads in turn to reduced emission line equivalent widths. 

4) A single, fixed, IMF is considered only by both groups. In view of the un- 
certainties on this quantity it appears useful to explore a wider range of 
IMFs. 

We here present results from new calculations which relax all the above as- 
sumptions and allow the exploration of a wide parameter space in terms of stellar 
tracks (including also alternate tracks with strong mass loss), IMF, and star for- 
mation history (instantaneous bursts or constant SFR) . A more detailed account 
is given in Schaerer (2001). The models are then used to simulate observations 
with a medium resolution cryogenic IR spectrograph for the VLT. 

2 Input physics and results for individual stars 

To account for non-LTE effects, which are crucial for a proper description of 
the spectra of hot stars considered here we use the TL US TY code of Hubeny & 
Lanz (1995). Several additional models were calculated with the CMFGEN code 
of Hillier & Miller (1998) to explore the possible importance of mass loss (as well 
known for Pop I stars, Schaerer & de Koter 1997) on the ionising spectra of Pop 
III stars. For the hottest stars (T c g > 80 kK) dominating especially the ionising 
spectrum in the He n continuum, no significant differences with plane parallel 
models are found as their continua are formed deep in the quasi static part of 
the atmosphere. 

We use stellar tracks from 1 to 500 M© calculated with the Geneva stellar 
evolution code (Desjacques 2000) or tracks from Marigo et al. (2001). Although 
probably unrealistic, we also explore the impact of high mass loss tracks from 
Klapp (1983) and El Eid et al. (1983). The atmosphere models and tracks are 
included in the synthesis code of Schaerer & Vacca (1998), which also calculates 
nebular line and continuous emission (adapted here to include additional lines 
and for an appropriate electron temperature T e = 30000 K). 

Based on these models we have calculated the ionising properties of ZAMS 
stars and their time averaged properties taking into account their evolution. The 
quantities are tabulated and given as fit formulae of use for various calculations 
(see Schaerer 2001). 

3 Properties of Pop III "galaxies" 

Integrated spectra have been calculated for a variety of IMFs. SEDs of integrated 
zero metallicity stellar populations are shown in Fig. [l] for the case of a Salpeter 



Observing the first galaxies with the VLT 



Pop III: Salpeter IMF (1-500 MJ 



Pop III: Salpeter IMF (1-500 M a ) 




2000 4000 

X [Al 




2000 4000 

x [A] 



Fig. 1. Spectral energy distribution (SED) including H and He recombination lines 
for Salpeter IMF from 1 to 500 M Q . Left panel: ZAMS population. The pure stellar 
continuum (neglecting nebular emission) is shown by the dashed line. For comparison 
the SED of the Z - 1/50Z Q population (Salpeter IMF from 1 - 150 M ) is shown by 
the dotted line. The vertical dashed lines indicate the ionisation potentials of H, He , 
and He + . Note the presence of the unique He II features (shown as thick dashed lines) 
and the importance of nebular continuous emission. Right panel: temporal evolution 
of the spectrum for 0, 2 and 4 Myr showing the rapid change of the emission line 
spectrum, characterised by the disappearance of the He II lines 



IMF from 1 to 500 M and instantaneous bursts of ages (ZAMS), 2, and 4 Myr. 
Overplotted on the continuum (including stellar + nebular emission: solid lines) 
are the strongest emission lines for illustration purpose. The striking differences 
when compared with non-zero metallicity, most importantly in the ionising flux 
above the He n edge (> 54 eV), have already been discussed by TS00. The 
comparison of the total spectrum (solid line) with the pure stellar emission 
(dashed) illustrates the importance of nebular continuous emission neglected in 
earlier studies (TS00, BKL01), which dominates the ZAMS spectrum at A > 1400 
A. The nebular contribution, whose emission is proportional to Q(H), depends 
rather strongly on the age, IMF, and star formation history. For the parameter 
space explored (see Schaerer 2001), we find that nebular continuous emission 
is not negligible for bursts with ages < 2 Myr and for constant star formation 
models. 

The spectra in Fig. [l] show in addition to the H and He I recombination lines 
the presence of the strong He n AA 1640, 3203, and 4686 recombination lines, 
which — due to the exceptional hardness of the ionising spectrum — represent a 
unique feature of Pop III starbursts compared to metal enriched populations (cf. 
TGS01, Oh et al. 2001, BKL01). Another effect highlighted by this Figure is the 
rapid temporal evolution of the recombination line spectrum. Indeed, already 
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Burst model, 10 M_solar, log(Age[yr])=4.0 



Burst model, 10 M_solar. log(Age[yr])=4.0 




Fig. 2. J (left) and K' (right) band magnitude as a function of redshift over the interval 
z ~ 6 to 14.5. Pure stellar continuum: upper dashed-dotted curve; intermediate/green 
dashed curve: including the Lyman-a, He II A1640, and He II A3203 emission lines; total 
predicted spectrum including lines and nebular continuous emission (lower curve/red). 
Note the increase by ~ 1 - 2.5 mag depending on z\ 



> 3 Myr after the burst, the high excitation lines are absent, for the reasons 
discussed before. 

The detailed behaviour of the relative line intensities and equivalent widths 
and their dependence on the IMF and the star formation history are discussed 
in Schaerer (2001). For the case of constant star formation, the He II or H 
recombination lines can be used as a measure of the star formation rate (SFR). 
For obvious reasons, strong variations of the various SFR indicators on the poorly 
known IMF are obtained. 



4 Simulations for a 2nd generation multi-object IR 
spectrograph on the VLT 

We have used the above synthetic spectra to calculate the photometric properties 
and to simulate spectroscopic observations. The following assumptions are made: 
A popular cosmology (i? m = 0.3, Q& — 0.7, H = 75 km/s/Mpc _1 ) is adopted. 
Lyman-a forest blanketing at redshift z < 6 is included following the prescription 
of Madau (1995). Sources are unresolved on a 0.3" scale, and the seeing adopted 
is 0.8" . We take typical SOFI and ISAAC near-IR filter responses to calculate 
magnitudes in the Vega system. Telescope parameters correspond to the VLT. 
The characteristics of the near-IR spectrograph are set similar to EMIR (Balcells 
et al. 2000), with 0.3" /pixel, 1" slit-width, and a mean total efficiency of 40%. 
All simulations shown here are calculated for a ZAMS population of 10 7 M Q 
with a Salpeter IMF from 1 to 500 M Q (somewhat more "favourable" than a 
constant star formation case), a total exposure time of 10 5 sec, and resolutions 
R ~ 1000 to 5000. 
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LYMAN ALPHA. Burst model, 10 7 ILsolar, log(Age[yr])=4.0 Hell 1640, Burst model. 10 7 M_solar, log(Age[yr])=4.0 




Fig. 3. Predicted S/N in Lyman-a (left) and He II A1640 (right) as a function z for a 
spectral resolution R — 1000 and t oxp = 10 5 s, for a 10 7 M Q burst at age ~ with a 
Salpeter IMF from 1 to 500 M 



As shown in Fig. [| the predicted magnitude in the K' band is ~ 24 - 25 
for a 10 7 M Q burst. This Fig. shows also the strong contribution of the nebular 
continuum, and the signatures of the strong emission lines on this broad band 
filter, which both need to be taken into account. Similar effects are observed in 
the other filter bands, the most relevant emission lines being Lyman-a and He n 
A1640. Simulations for various cases of the IMF show that standard broad-band 
colors do not allow a distinction between upper mass cut-offs of 100 to 1000 M© 
and are also insensitive to changes of Mi ow . 

The expected S/N in Lyman-a and He II A1640 as a function of redshift 
of the source (observed in the JHK bands) for a spectral resolution of 1000 is 
plotted in Fig. H. The calculations illustrate the following: 

• Lyman-a can easily be detected with a good S/N over a redshift ranges (z ~ 
8-10.5, 11.5-14, and 15.5-18). A joint detection with He II A1640 is possible 
for z ~ 5.5-7.5 (Lyman-a in optical domain) and z ~ 8-10; at larger redshift 
the S/N of He n A1640 becomes very low. 

A detection of both He n A1640 and Lyman-a is obviously important to 
secure the redshift, and to obtain a measure of the hardness of the ionising 
flux. This allows potentially the distinction between different IMFs. 

• Higher spectral resolution (R ~ 5000) considerably increases the chances of 
detection between the sky lines. For unresolved lines the resulting decrease 
of S/N is modest. The medium spectral resolution is also required to attempt 
to measure the emission line profiles, in order to distinguish Pop III source 
from AGN (cf. TGS01). Once this is obtained, the data can be rebinned to 
increase the S/N. 

How many such Pop III galaxies are expected ? The number of metal-free 
objects depends of course on the galaxy formation scenario, the star formation 
efficiency, and the hydrodynamical enrichment process of the ISM, among other 
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factors. The comprehensive study of Ciardi et al. (2000) shows that at z > 8 
naked stellar clusters, i.e. objects which have complete blown out their ISM 
dominate the population of luminous objects. One may expect that all metals 
are ejected in this way into the IGM, thereby avoiding largely local pollution, i.e. 
leaving these objects as Pop III "galaxies". Based on such an assumption and 
using the ionising fluxes from TS00, Oh et al. (2001) have recently calculated 
the predicted number of Pop III starbursts detectable in He n lines with NGST 
for a one day integration. Their estimate yields between ~ 40 and 3000 sources 
in a 8'x8' field of view envisaged for the proposed IR multi object spectrograph. 
Although pilot studies have recently started to explore the formation of dust in 
Pop III objects (Todini & Ferrara 2001), the effect is neglected here. 

In conclusion, a medium resolution (R up to ~ 5000) IR multi-object spec- 
trograph on the VLT covering the J, H, and K band would not only be of 
great interest to numerous studies of intermediate and high-redshift galaxies (cf. 
contributions of Cuby, Genzel, Lefevre, this meeting). It should also allow obser- 
vations of the first so called Pop III objects — certainly one the great remaining 
challenges of observational astronomy. If realised on a reasonable time scale, such 
ground-based observations with the VLT at ESO should be able to achieve an 
important breakthrough before the launch of NGST. 
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